Backward uniqueness for the heat equation in cones 
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Abstract 

It was shown in [5l [13] that a bounded solution of the heat equation 
in a half-space which becomes zero at some time must be identically 
zero, even though no assumptions are made on the boundary values of 
the solutions. In a recent example, Luis Escauriaza showed that this 
statement fails if the half-space is replaced by cones with opening angle 
smaller than 90°. Here we show the result remains true for cones with 
opening angle larger than 110°. 

1 Introduction 

Consider an open set C M". Let u be a bounded solution of the equation 

ut- Au + b{x,t)Vu + c{x,t)u = inilx(0,T), (1) 

where the coefficients 6 = . . . , 5„), c are measurable and bounded. We 
say that has the backward uniqueness property if the following statement 
holds: 

(BU) If a bounded u: 17 x (0, T) ^ M satisfies © and u{-,T) = 0, tlien 
u = inQx (0,r). 

It is important to emphasize that no assumptions are made about u 
at the parabolic boundary d^l x (0,T) U (0 x {0}). In fact, we can think 
about the problem in terms of the control theory: we are given some initial 
data tto : 17 — 5- M, and we wish to find a suitable boundary condition g on 
the lateral boundary di} x (0, T) so that when we solve equation ([T]) with 
uq as the initial condition and g as the boundary condition, the solution 
will become exactly zero at time t = T. In this interpretation condition 
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(BU) means that we can never achieve the exact boundary control of any 
non-trivial solution. 

While the control theory for PDEs seems to be the most natural back- 
ground for (BU), the problem also appeared in regularity theory of parabolic 
equations, such as the Navier-Stokes equations, harmonic map heat flows, or 
semi-linear heat equations, see [5l[TTl[l3]. The specific unbounded domains 
which arise in this connection are complements of closed balls (for interior 
regularity), half-spaces (for boundary regularity at boundaries), or cones 
(for boundary regularity in Lipschitz domains). 

By classical results we know that in bounded domains we can achieve ex- 
act control, and therefore any domain satisfying (BU) has to be unbounded. 
Classical backward uniqueness results for parabolic equations imply that 
Q = satisfies (BU). It turns out the the half-space = R" also satisfies 
(BU), although this is harder to prove, see [5]. In general, the smaller the 
domain, the harder it will be to show that it satisfies (BU). It is immediate 
that if fii C and fii satisfies (BU), then also satisfies (BU). 

In this paper we consider the question for cones with opening angle 6. 
In suitable coordinates 

Og = {x = {xi,x'),x' £ M^-^xi > |x|cos(6l/2)}. (2) 

Luis Escauriaza [1] recently showed that - surprisingly - (BU) fails when 
9 < 7r/2. We briefiy recall the counterexample. Let us denote by T the stan- 
dard heat kernel, i. e. T{x,t) = (47rt)~"/^ exp (|2;p/4i), and recall Appell's 
transformation 

X 1 

u{x,t) = T{x,t)v{y,s), y = -, s = j. (3) 

This transformation takes the solutions u{x, t) of the heat equation into 
the solutions v{y,s) of the backward heat equation 



Vs + Av = . (4) 

By taking u{x,t) = h{x) for a suitable harmonic function h in Og, we can 
get a counterexample to the backward heat equation form of (BU) for 6 < 
tt/2. In dimension 2 we can use the real part of the holomorphic function 
z —7- exp{—Az") (for suitable ^ > and a parameter a > 2) to obtain an 
explicit formula: 

= Re iexp(-A<^^ll^ + !#). (5) 



2 



The function v{yi,y2, s) is bounded in a sector that | arctan < 7r/(2a), 
away from the origin. We can shift it to f (yi + 1, 2/2) s). The resulted function 
is bounded in a sector with angle n/a satisfying the backward heat equation 
Vs + Av = 0, and v{-, -,0) = 0. 

We note that it is enough to construct a counterexample in dimension 
n = 2. The higher-dimensional example can then be constructed by sim- 
ply considering the two-dimensional function as a function of n variables, 
independent of X3, . . . , x^. 

Escauriaza's example shows that (BU) fails for 9 < 7r/2. Since we also 
know that (BU) is true for = vr, it is easy to see that there exists a 
borderline angle £ ['^/2,7r] such that (BU) is true for 6 > 9q, and (BU) 
fails for 9 < 9q. The borderline case 9 = 9q might perhaps present an extra 
difficulty. 

The main result of this paper is the following: 
Theorem 1.1. The cones Og satisfy (BU) for 

9 > 2arccos(l/\/3) ~ 109.52° . 
In other words, the critical angle 9q introduced above satisfies 

9o < 2arccos(l/\/3). 

It is tempting to conjecture that 6*0 = vr/2. This is supported by the fact 
that = 7r/2 is the borderline case for the above construction of Escauriaza, 
as can be seen from the Phragmen-Lindolef principle. 

For the classical heat equation, corresponding to the case 6 = and c = 
in ([1]), and 9 = n (the half-space), the statement (BU) can be proved by 
a relatively simple application of Fourier transformation and some classical 
complex analysis results, see [13]. We were not able to find such simple proof 
of the case b = 0, c = when 9 < tt. 

While completing this paper, we learned about reference [8jll] Theorem 
6 in [8] states that for the classical heat equation (corresponding to 6 = 
0, c = 0), (BU) holds if and only ii 9 < tt/2. Unfortunately, it seems the 
proof is not available in print. 

Our proof of Theorem (jl.ip relies on two Carleman-type inequalities, 
along lines similar to [5]. The first inequality. Proposition 12. 11 is taken from 
[5] and is applied in the same way to obtain fast decay rates for the solutions, 
see Lemma l2.2i We note that Carleman inequalities of this type are can be 
found already in [21 El El Sj . 

^We thank Gregory Seregin for pointing out this article. 
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The second inequality, Proposition 12.3^ is the main new tool used in 
our proof. The heuristics behind this inequality is somewhat similar to the 
heuristics behind the second Carleman-type inequality in [5] (Proposition 
6.2). However, the proof of Proposition 12.31 requires a new idea, as for 
6 < TT certain critical terms appearing in the proofs of the inequalities lose 
convexity. 

In addition to determining the critical angle, another interesting open 
problem is to optimize the assumptions on the coefficient b{x,t) and c{x,t), 
in the spirit of [9]. For example, it is conceivable that the result remains 
true for b G L'^'t'^ and c G L^"^^^''^, but it might be a difficult problem to 
decide whether this is the case. 

In what follows we will work with the inequality 

\ut - Au| < ci(|Vm| + (6) 

rather than It is not hard to see that when assuming the boundedness 
of b and c, the two formulations are equivalent. 



2 Backward uniqueness 

Without loss of generality we assume T = 1 and work with the backward 
form of dll). 
Recall that 

Oe = {x = {xi,x'),x' G M"-\xi > \x\cos{9/2)}. (7) 

Suppose that u{x,t) is a solution to the backward heat equation for some 
e > 2arccos(l/A/3). 

\ut + Au\ < ci{\Vu\ + \u\) inC'0x(0, 1), (8) 
u{-,0) = in Oe. (9) 

In addition, 

|m| < M in Oe x(0, 1). (10) 

Then u = 0. 

To prove the above statement we firstly need the following Carleman in- 
equality from [5] , by which we obtain a decay result for solutions of backward 
heat equation. 
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Proposition 2.1 ([5J). For any function u G C^(M"' x (0, 2);]R") and any 
positive number a, 

[ /i-2«(t)e-^ (-|n|2 + \Vu\^) dxdt 

<co [ /i"2"(t)e"^|9jti + Anl^dxdt, (11) 

JR"x(0,2) 

where cq is a positive absolute constant and h(t) = te~ . 

Lemma 12.21 below immediately implies exponential decay of the solution 
u. The proof is by using the Carleman inequality in Proposition 12.11 The 
decay of u enables us to apply the Carleman inequality in Proposition 12.31 
and reach the conclusion in Theorem II. 1[ 

Lemma 2.2. Let Bn denote the ball with radius R in M". Assume that 
R > 2. Consider a function u satisfying the following conditions, with some 
positive constants ci and M. 

\ut + Au\<ci{\Vu\ + \u\) inBRx{0,T), (12) 

n(x,0) = in Br, (13) 

|n| < M in Br x (0,r). (14) 

Then there exist some constants (3, 7, such that for t G (0,7), 

um< . f ^. Me-^'r, (15) 
mmj 1 , i I 

where j3 is a small enough absolute constant, C2 depends on ci, 7 depends 
on ci and T . 

We will give the proof of the lemma in the next section. 
The following Carleman inequality in Proposition 12. II is a key tool used 
in our proof of the backward uniqueness in cones. We define the set 

Qe = {Oer\{xi > 1}) x (0,1). 

The purpose of "cutting off the corner" is to avoid singularities at the origin. 

1 — t 

Proposition 2.3. Let 4>{x,t) = aK{t)ip{x) + t^ , where A(t) = ^^^^ , and 

(p{x) = xf — 6°' r°, where r = \x\, and e = cos{6/2). For any e G (0, l/\/3), 
that is, G (2 arccos(l/\/3), vr), there exists some a = a{e) G (1,2) such 
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that the following inequality holds for u £ C^{Qg) and a > oq for some 
constant oq. 

[ e^-^^^'*) [a (A(t) + ifix)) + | Vn|2] dxdt 

<a[ e^'^^'^'^^ldtu + Aufdxdt. (16) 

We apply this Carleman inequality to prove the main result of the pa- 
per in the remaining part of this section. The proof of Proposition 12.31 is 
postponed to the last section. 

For X £ Og we denote by dg^x) the distance between x and the boundary 
of Og, explicitly given by 

dg{x) = xi sm{e/2) - \x'\ cos(6l/2). (17) 

Let Og^'^ = {x e Og \ dg{x) > 2}. With any other number c, the set Og'^'' 
is defined in the same way. 

The next lemma is a consequence of the decay result from Lemma 12.21 
and Proposition 12.31 . It implies Theorem 11.11 immediately. 

Lemma 2.4. Assume that for some 9 G (2 arccos(l/-v/3), ti") a function u 
satisfies (0) - J70|). then there is a number 71 (ci) G (0,7/2) such that 

u{x,t) = Q (18) 

in x(0,7i). 

Proof. Lemma 12.21 implies that 

\u{x,t)\<C2Me-l^— (19) 

for all {x,t) G Oe"*"^ x(0,7). By local gradient estimates for the heat equa- 
tion [ IQj we can assume that 

|n(x, t)\ + \Vu{x, t)\ < csMe 2?— (20) 

for all {x,t) G 0^+3 x (0, 7/2]. 

By scaling we define a function v by 

viy,s) = u{Xy,Xs^ -^i) (21) 
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for (y, s) & O0 X (0, 1) with A = -^271. This function satisfies the relations 

\dsv + Av\ < ciX{\Vv\ + \v\) inC'0x(O,l) (22) 
viy,s) = in 0,x (0,1/2), (23) 



and 

|v(2/,s)| + |Vz;(y,s)| < csMe 2(a^--7i) < csMe"^^ (24) 

for 1/2 < s < 1 and y G 0^+^/^ = {y £ Og \ de{y) > 3/A}. 

To apply Proposition 12.31 we need certain decay of \v{y,s)\ when |y| 
is large. Notice that the preferred decay can be obtained by considering 
a cone with slightly small opening. Proposition 12.31 holds for angles in 
(2 arccos(l/-v/3), vr). We thus consider the median of 6 and 2 arccos(l/\/3), 

, 6' + 2arccos(l/\/3) 

= . 

2 

In the smaller cone Os = {x £ > |x| cos((5/2)} we have the estimate 

deiy) > \y\ sin (^). It follows ([Ml) that 

\v{y, s)\ + \Vv{y, s)\ < csMe"^'^ (25) 

for 1/2 < s < 1 and y G O^nOe+^Z^, with the constant /3' = /3sin2(^). 
We can further have 

\v{y, s)\ + \Vv{y, s)\ < ^Me"^'^ (26) 

for 1/2 < s < 1 and y £ Os (^{yi > 3/A} for some other constant C3. 

Next we work on the smaller cone Os with opening 5, where we have 
exponential decay (|26|) and the following properties inherited from Oq. 

|(9^u + A?;| < ciA(|Vt;| + |f|) inC'5x(0,l), (27) 
v{y,s)=0 in O5X (0,1/2). (28) 

Proposition 12.31 requires support condition for the Carleman inequality. 
For that purpose, let us fix two smooth cut-off functions such that 



V'i(yi) = 

V'2(t) 



0, yi < 3/A + 2, 

1, yi > 3/A + 3, 

0, T < -3/4, 

1, T>-l/2. 
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We set (for the definition of (f), see Proposition I2.3P 

My, s) = -Hy, s)-B = {i- + '--B, 

a s"'^ a 

where e = cos{5/2), B = ^4>{yx, with yx = (3/A + 3, 0, • • • ,0) and 

V{y, s) = ^i(?/i)V'2(-^), w{y, s) = ??(y, s)v{y, s). 

The function w is not compact supported in Qs (recall that Qs = {Os n{xi > 
1}) X (0, 1)). However, it follows from (|24p and the special structure of 
the weight function (p in Proposition 12.31 that, with w replacing u in (|2.3p . 
integrals on both sides converge. If we multiply w by an additional cut-off 
function ^ such that 

1, T <R 
0, T>2R 



e(x) 



and |V.^| < c/i?, |V^.^| < c/B?, apply Proposition 12.31 to the compact sup- 
ported function w^, then let R — )• c«, we finally obtain 

f e^'"^^ [a (A(s) + ^)w^ + \Vw\^] dyds 

< 4 /" e^'"l'B\dsW + /^w\^dyds. (29) 

From (12711 we have 



\daW + /\w\ < ciA(|Vt(;| -I- 

+C4(|Vz;| + \v\){\dsV + |V7?| + |A77|). (30) 

Refer to the definition of ip2 we know that (ps/B > —3/4 in the support of 
w. In other words aA{s)(p{y) -|- > aB/A. For a large enough, this implies 
that a(A(s) -|- (p{y)) > 1 by Cauchy-Schwartz inequality. In addition, we 
take 7i(ci) small enough such that WcfX^ < 1/2. We then have 

I = [ e^""^^ {w"^ + \Vw\^])dyds (31) 
JQs 

< 32cl I e^'"f'^{\v\^ + \Vv\^){\d,v\ + \Vr]\ + \Ari\fdyds. (32) 

To estimate the right hand side, we need to look into the detail of deriva- 
tives of r/. In view of the definitions of ipi and ■02; the support of derivatives 
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of ri{y, s) is the closure of the set 



{yi > 3/A + 2, -3S/4 < 0b < -B/2} 

U {3/A + 2 < yi < 3/A + 3, > -B/2}. 

However, the second set has empty intersection with Os x(l/2, 1), where 
the function v is nonzero. Hence the support of the term (IV^I + + 
iVr/l + |A?7|) is closure of the set 

a; = {yi > 3/A + 2, 1/2 < s < 1, -SB/A < (j)B{y, s) < -B/2}. 

We denote by xiu^ the characteristic function of the set co. 

Next we estimate the rate of increasing at infinity of derivatives of r]{y, s) 

1 — s 

in the set uj. Recall that (/>(y, s) = aA{s)ip{y) + s^, where A(s) = — — - and 

ga/2 

viv) = Ui The function A(s) and the derivative A'(s) are bounded 

for s G (1/2, 1). The function if and its derivatives up to the second order are 
bounded by (const. The cut-off functions ipi and '02 and derivatives 

up to the second order are bounded by some absolute constant. The value 
of B is bounded from below regardless of the value of the parameter a. Thus 

(|9,r?| + |Vr?| + |Ar/|)2<c5|y|2° (33) 

in the set uj. 

We now estimate ([32]) by using ([25]) and ([33]) . We see that 

/<C6Me-^» / |y|2-e-2^'^x(y,5)dyd5 
■JQs 

for some constant cq. The last integral in bounded. Passing to the limit as 
o — )• DO we see that v{y, s) = for 1/2 < s < 1 and 0B(y, s) > 0. Using the 
property of unique continuation across the spatial boundaries (see Theorem 
4.1 in we show that v{y,s) = if y G and < s < 1. This proves 
the lemma. □ 

3 Proof of Lemma 2.2 

The proof of Lemma 12.21 is based on the Carleman inequality in [5] , which 
we quoted in Proposition 12. li 

Proof of Lemma \2.S\ The proof is similar to the one in [5] . We still include 
the proof here for the convenience of the reader. 
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In what follows, we always assume that the function u is extended by 
zero to negative values of t. 

The assumption that R > 2 results in no loss, since the conclusion is 
only useful when R is large. According to the local gradient estimates of the 
heat equation [TU], in the smaller cylinder {x,t) E x (0, r/2), we can 

assume that 

\u\ + \Vu{x,t)\ < . M (34) 

mmjl, 1 ) 

with some absolute constant cy. 

We fix t G (0, min{l, T}/12) and introduce a new function v by the usual 
parabolic scaling: 

v{y,s) = u{Xy,X^s -t/2). 

The function v is well defined on the set Qp = B{p) x (0,2), where 
/9 = (i? - 1)/A and A = \/3t G (0, min{l, \/T}/2). We have the following 
relations for v. 

I^sf + A^l < ciAdVi-l + 1^1), (35) 

\v{y, s)\ + \Vv{y, s)\ < —^L—M (36) 

mmj 1 , i j 

for all {y,s) G Qp, 

v{y,s) = ^ (37) 

for {y,s)eB{p) x (0,1/6]. 

By the assumption that R > 2 and A < 1/2, we have p > 2. In order to 
apply Proposition 12. H we take two smooth cut-off functions in Qp-. 

( 0, \y\ >p-l/2, 

^^^^^ = 1 1, \y\<P-i, 



Ms) 



0, 7/4 < s < 2, 

1, 0<s<3/2. 



By assumption, these functions take values in [0, 1] and are such that | V^^/^pl < 
Ck, k = l,2, and {dgiptl < Cq. We set ri{y,s) = ipp{y)'il^tis) and 

w{y^s) = r]{y,s)v{y,s). (38) 

It follows from ()35p that 

\dsw + Aw\ < ciA(|Vu;| + \w\) + csx{\^v\ + \v\), (39) 

where cg is a positive constant depending only on ci and C^, k = 0,1,2, 
Xiy,s) = 1 for {y,s) £u; = {p-1 < \y\ < p,0 < s < 2}U{|y| < p- 1,3/2 < 
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s < 2}, and x{y-,s) = for (y, s) ^ uj. The set oj is where the cut-off 
function r/ is not constantly 1 in Qp. Obviously, the function w is compactly 
supported on x (0,2), we may apply Proposition 12.11 and obtain 



/ /i-2-(s)e-^f-|H' + |V^ 

J On 



w\ ) dyds 



< 



Co 



h-'"'{s)e-^\dsw + /^wydyds. (40) 



Taking a > 2, and applying (I39p we finally obtain that 



where 



Qp 



h-"'{s)e-'^{\w\^ + \Vw\')dyds <Aco{ciyi + cih), (41) 



= [ x{y,s)h-'^''{s)e--^{\Vv\'^ + \v\'^)dyds. 

Jqp 



Taking a sufficiently small value for 7 = 7(ci) such that in the range A £ 
(0,7), we can assume that the inequality Acqc\\^ < 1/2 holds, and therefore 
(jlT]) implies that 

/ < Sc^clh. (42) 
Notice that near the origin {y = 0, s = 0}, where the parametric function 

|y|2 

/i~^"(s)e~^ is not integrable, our characteristic function x is 0. By ([36]) 
we have 



h < 



min{l, r^} 

f2 



h ^"(s)e dyds 



3/2 J|y|<p-i 



< 



Jp-l<\y\<p 
C9M2 



/i (s)e 4s dyds 



min{l,T2} 
where cg is an absolute constant. 



/i-'''(3/2)+ / /i-"''(s)e-"-4^ds 



(43) 
(44) 
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Using ()44p we obtain the estimate 



D = 



B{1) Jl/2 



\w\ dyds 



B{1) Jl/2 



\v\ dyds 



< cio / h'^''{s)e~'^{\w\'^ + \Vw\'^)dyds 

hp 



< 



cuM^ 



min{l, T^} 

cuM^ 



/i-'"(3/2)+ / h-"'{s)e-^ds 



min{l, T^} 
We take (3 < 1/64 and then let 

a = /3pV(21ogM3/2)). 
This choice of a leads to the estimate 



D < ciie-^p" 



1+ / g{s)ds 
Jo 



where g{s) = h {s)e 32s. By simple calculation we have that 



g'{s) = h-'''{s)e-^ 



log/i(3/2)^s 3^ 32s2 



One can readily verify that g{2) < 1 and g'{s) > for any s € (0,2) if 
/? < g^log/i(3/2). Therefore, 



D < 3 -^—^e 



'^p' = 3- 



min{l,r2} min{l,T2} 
On the other hand, the regularity theory implies that 
|M(0,t)|2 = |t;(0,l/2)|2 <ci2D. 

Finally we obtain that 



e ^ i2t . 



u 



{0,t)\<—^—-Me-^^^. 
mmj 1 , i I 



Taking another constant equals /3/24, and still denotes by /3, we reach the 
conclusion of the lemma. □ 
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4 Proof of the Carleman inequality 

One of the difficulties in the proof of the Carleman inequality in Proposi- 
tion 12.31 is that - by comparison with the case 9 = n - some loss of the 
convexity of the weight ip cannot be avoided. Therefore we have to investi- 
gate in more detail some of the terms which can be neglected when Bq > ir. 

Proof of Proposition We denote (f) = aA{t)(p{x) + t'^. Let u be an arbi- 
trary function in C^{Qe) and v = e'^u. Then 

Lv = e'I'idtu + Au) = Av + |V(/>|^v - dt4>v + dtv - 2V(j)Vv - A<j)v. (45) 

We decompose L into symmetric and skew-symmetric parts 

L = S + A, 

where 

Sv = Av + \V^\'^v-dt(l)v (46) 

and 

Av = dtv - 2V(j)Vv - A(j)v. (47) 
The right hand side of the inequality (jl6p is 



Lv\'^dxdt = J \Sv\'^dxdt + J \Av\'^dxdt + J {[S,A]v)vdxdt, (48) 

where [S, A] = SA — AS is the commutator of S and A. By simple calcula- 
tions we have that 

{[S,A]v,v) = j A(j)^kiv,kV,idxdt (49) 

+ j (2V(/.V|V(/)p - A'^cl) + d'^cp - 2dt\V(l)\'^) \v\'^dxdt.{50) 

The Hessian of the function cp = aA{t){x'^ — e" r") is not positive-definite. 
To compensate the term J 4:(j)^kiv,kV,idxdt in (liOl) we introduce a function 
F{x, t) to be determined. 

{Sv, Fv) = J AvFv + (I V(/>|2 - dt(t>)Fv'^dx 

= j -FjVup + i^AF + |V(/>|2f - dt(t>F)v^dxdt. 
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Cauchy-Schwartz inequality implies that 

{Sv,Sv) > -{Sv.Fv) j F'^v'^dxdt 

> j F| Vup - AF + \V<t)\^F - dt(t)F + ^F'^)v'^dxdt. (51) 
Combining ()49p and (I5ip we have 



{[S, A]v, v) + {Sv, Sv) > J 4:(j),kiv,kV,i + F\Vv\^dxdt 
+ j {iVcfiVlV^l"^ - A'^(t> + d^(t>-2dt\V(t>\'^)v'^dxdt 
+ j -{^AF + 1^4)1"^ F-dt(t>F + ^F^)v^ dxdt. 
By calculation the Hessian of (f is 



(52) 
(53) 
(54) 





f xr' 


•• 


• M 







•• 


• 


D^(p{x) = a{a — 1) 










V 


•• 


• 0/ 



ae"r"-2F„ (55) 



+a(2-a)e" r'^'^x^ x, 



(56) 



where En reprensent the n dimensional identity matrix, x = {xi, . . . ,x„) is 
the row vector and x'^ denotes the transpose of x. It is easy to see that 

D^^{x) +ae''r'^-^En > 0. 

We thus let f{x) = ae" r"~^ and 

F{x,t) = AaA{t)f{x) + l. (57) 

With this choice of F(x,t), the right hand side of line (]52p is positive and 

4>,kiv,kV,i + F\'Vv\'^dxdt > J \'Vv\'^dxdt. 

Grouping the remaining terms according to the orders of the parameter a, 
with A3 denoting the terms with and etc, we have that 



{[S,A]v,v) + {Sv,Sv) > 



iVvl dxdt 



+ j {A3 + A2 + Ai+ Ao)v^dxdt (58) 
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where 



A3 = 4a^A-V,fci95,fcV',« - 4a^AV|Vv?|^ (59) 
A2 = -4a2AA'|V9?p + 4a2AAV/-4a2A2(t)/2-a2A2|V(^P, (60) 
Ai = -aAAV + aA"(p - 2aA(t)A/ + aAV - 2aA(t)/ + 8atAf, (61) 



Ao = 7/4 + 2t. 



(62) 



We analyze A3 first. With as a coefficient A3 must be non-negative in 
the set Qe. By letting xi/r ^ e, we see easily that e < l/-\/3 is a necessary 
condition for A3 > 0. This could be seen by letting . Next we show that 
under the condition e < l/VS, we indeed have that A3 > 0. Denoting by 
Vip'^ the transpose of the row vector Vi^, we notice that 



A3 = 4a^A^{t)V^ (-DV(^) - a£"r"-2) V<^^. 



(63) 



It is easy to see that 

( 



„a-2 












> a{a — 1) 

V 

Apply the fact that > r"^~^, we have 



0\ 




2q;£" r"-2 



-En 



(64) 



a ■ 



2£° 









-2 £«£;„_! 



(65) 



where E^-x is the n — 1 dimensional identity matrix. 
The first derivatives of ^ are as follows. 



Xk, 



k = 2, . . . ,n. 



We notice that 99,1 > a{l - e°')x'^~\ Thus 
A3 > 4o3A3(i)aV"-2 [{a - 1 - 2£")(1 - e 



a\2 2o-2 



X 



2g3a^2a-4, /|21 



Taking into account that xi/r > £ and |a;'| /r < 1 — £ , 

A3 > 4a^A3(t)aV"-^£2"-2 [(a - 1 - 2£°)(1 - e^'f - 2e''+'^{l - e'^)] 
Let us denote the quantity in the bracket above by m{a,e). 



m{a,s) = (a - 1 - 2£")(1 - £«)" - 2£"+2(l - s') 
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^3 is non-negative if m{a, e) is. In the set a G (1, 2) and e G (0, l/\/3), the 
function m{a, e) is monotone increasing with respect to a and is monotone 
decreasing with respect to e. Notice that 



m 



(2,l/\/3) = 0. 



Hence for any e < l/\/3, there exists a corresponding a{e) < 2 such that 
for a £ (a(e), 2), m{a, e) > 0, and in turn A3 > 0. 

For the estimate of A2 we notice that Vv = V0f + e'^Vu, a.s v = e'^u. To 
bound |e'^Vn|, we want to apply the inequahty |e'^Vnp/2 < |Vt'p + |V(/)pf ^. 



We thus look at the inequality (|58p in the following way. 

{[S,A]v,v) + {Sv,Sv) > [{\Vv\'^ + \V^\'^v^)dxdt 



(66) 



+ j[A3 + {A2 - |V0|2) + Ai + Ao]v^dxdt 

Next we estimate A2 — |V0p. 

A2 - ^ -,a^Aimt) ((1 + ^) iv.i^ - ./ + 

A{t) = and A'{t) = . "Z^+jW^^ |A(t)/A'(i)| < ^. 

^2 - |V,^|2 > -4a'A{t)A'{t) ((1 - ^) |V^|' - 9^/ - ^/') . (67) 

Notice that jVc^p > |^^i|2 > ^^(i _ £a)22,2a-2 ^ ae"(l-e")x^"-2. 

Prom the expression of m{a, e) in the estimation of A3 above, we know that 
e < (a — l)/2. Taking into account that r > xi > 1, 

A2 - \V(t>\^ > -Aa^A{t)A'{t) Qx?""' - ^^'""'') > -yA(t)A'(t)x2°-2. 

(68) 

Finally, we estimate Ai. Recall that 

Ai = -aAAV + aA"ip - 2aA(t)Af + aA'ip - 2aA{t)f + 8atAf. 
An simple observation is that 

Ai > a{A"ip + aA')ip - aA(AV + 2A/ + 2/), (69) 
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and A"(i) + A'{t) > a — 1 > 0. The terms in the second parenthesis are of 
homogeneity less than a — 2, thus under the control of A2. Consequently, 
there exists some constant oq depending on ip such that for a > uq, 



A2 + A^> -^A(i)A'(t)a;2"-2 + a{a - l)ip. 



(70) 



In addition |A'(t)| > 1 and xi > 1. We thus have 

{[S, A]v, v) + {Sv, Sv) > y (I Vf p + \V(p\'^v^)dxdt 



+ 



J '^K{t)v^dxdt + J a{a- l)ipv-' 



dxdt 



In turn 



/ 

JQ, 



^2aA{t)tp{x)+2t'^ 



— A(i) + a(a - 1)^ j + ^\yu\^ 



dxdt 



jQe 

To simplify, we can assume a > 3/2, and a > 2, it follows that 

/ e^-HtM-)+2t' [a (A + ^)u^ + \Vu\^] dxdt 
JQe 

< 4 / e2«^W^(^)+2*' \dtu + Aufdxdt. 
jQf, 



(71) 



(72) 



□ 
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